International Journal of Electrical and Computer Engineering (IJECE) 
Vol. 9, No. 1, February 2019, pp. 66~77 
ISSN: 2088-8708, DOI: 10.11591/ijece.v911.pp66-77 o 66 


Implementation of a frequency control in 
a biomass gasifier system 


Yecid Mufioz!, Adalberto Ospino?, Carlos Robles?, Carlos Arizmendi‘ 


^Facultad de Ingeniería, Universidad Autónoma de Bucaramanga, Colombia 
?Facultad de Ingeniería, Universidad de la Costa, Colombia 
Facultad de Ingeniería, Universidad del Magdalena, Colombia 


Article Info ABSTRACT 
Article history: Distributed power generation has grown in popularity in recent years, 
l especially in areas not connected to the national grid. As a result, rural 
Received Feb 6, 2018 microgrids are becoming more common, involving great potential for energy 
Revised Aug 10, 2018 based on biomass conversion such as gasification. After analyzing distributed 
Accepted Oct 8, 2018 power generation facilities in developing countries, the authors found 
—_——_— problems with the frequency stability. This paper focuses on solving the 
Keyword: problem of frequency control in energy supplied by microgrids based in 
biomass gasification. For that purpose, the authors have developed a physical 
Biomass model of a downdraft gasifier, this model was used for design a novel 
Downdraft Gasifier strategy for frequency control, which has been based and validated on an 
Microgrid existing gasification system, which supplies power to a population in Necocli 

(Colombia). 


Copyright © 2019 Institute of Advanced Engineering and Science. 
All rights reserved. 


Corresponding Author: 


Adalberto Ospino, 

Facultad de Ingenieria, 

Universidad de la Costa, 

Calle 58 No 55-66, Barranquilla, Atlantico, Colombia. 
Email: aospino8@cuc.edu.co 


1. INTRODUCTION 

The concept of distributed power generation (DG), or microgrids that work in isolation, is currently 
becoming an option to solve the problem of electricity supply, especially in rural areas in developing 
countries where the main grid is unavailable [1]-[4]. DG involves improvements in voltage profiles and 
reduces distribution losses because the generators are located near the rural load centers [5], [6]. 

The feasibility of biomass gasification in the rural electrification of remote areas has been widely 
investigated [7]. In [8] are considered different circumstances using an annual cost model taking into account 
the time-value of money, and suggests gasifier mass production as a prerequisite for making this technique 
feasible. Technology is currently closer to fulfilling this prerequisite and considering the rising price of fossil 
fuels and electricity, gasification is more applicable. Other authors present rural electrification as a favorable 
scenario for the gasification of biomass [9]. Biomass has a high potential contribution to the renewable 
energy mix in developing countries [10]-[12]. In the case of India, biomass-based power is widely available 
and more uniformly distributed than other renewable energy sources [13]. The potential for generation 
through renewable sources in India is 85GW and biomass could supply approximately 25% (about 20GW) 
and so represent a highly viable solution for the energy demands of small villages [14]. 

As in conventional grids, to maintain quality levels the frequency of the electricity must be 
maintained within a minimum deviation when using DG [15]. The authors in [16] suggest that DG units must 
contribute to frequency control and in order to contribute to primary frequency control these units must 
increase their power output. Various studies present simulations using the model of a conventional power 
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plant [17], [18], wind turbine [19], micro-turbine [20], fuel cells [21] and spark ignition engine power 
system [22]. 

The Universidad Politécnica de Valencia (UPV) and the Instituto de Planificación de Soluciones 
Energéticas de Colombia (IPSE) examined biomass generation systems in rural off-grid areas of Colombia. 
They found serious power problems in the frequency stability of the network, with variations greater than 
+10% of the desired frequency of 60 Hz. These problems are mostly caused by the lack of automation of the 
grid frequency control in biomass units. 

To solve this problem, this paper firstly proposes a physical model of the dynamic behavior of an 
isolated downdraft biomass gasification system which will be used to design a cascade control system. This 
model is based and validated on an existing gasification system in Necoclí (Colombia). A global analysis has 
been made. Results from other authors, experimental data, and new contributions have been used to develop 
the physical modelling. This modelling takes into account the nonlinearities of the analyzed systems. 


2. SYSTEM PHYSICAL DESCRIPTION 

The physical system consists of a three-phase synchronous generator driven by a reciprocating 
internal combustion engine (RICE) powered by synthesis gas provided by a downdraft biomass gasifier. 
Figure 1 and Table 1 shows images of the system and gasifier specifications respectively. 


Figure 1. RICE-generator (a) and (b) Gasifier installed in Necocli (Colombia) 


A schematic diagram of the system is shown in Figure 2, which also identifies the components. 
The biomass is stored in the hopper and fed to the gasifier through the double door feeder system. A limited 
amount of air for partial combustion enters through the air valve to the combustion zone to ensure the 
production of a relatively clean and good quality gas. 
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Figure 2. System layout 


The reactor maintains a coal bed that enables the reduction (or gasification) of the partially burnt 
products, and facilitates the evacuation of ashes through a perforated wall (pushed into the ash collection 
chamber when necessary). The gas outlet is connected to the power plant through a Venturi Scrubber, 
followed by coarse filters, fine filters, and a security filter. These filters are then followed by a security flame 
and the output valve of the plant (for gas cleaning and operation). 
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After the output of the gasifier, the synthesis gas (hereinafter, gas) generated by the gasification 
system, passes through a gas flow control valve (item 11 in Figure 2) and is sucked into the blower (item 12 
in Figure 2). Finally, the gas goes through the engine valve to access an internal combustion engine coupled 
to the generator. In the case of the gasifier under consideration, air enters through the air nozzle of the 
hopper, and the amount entering the generator depends on the blower speed and the gas flow control valve. 
Specifications of the blower and RICE-generator are presented in Table 2. 


Table 1. Gasifier Specifications 


Characteristic Value 

Model Ankur WBG-80 

Operating mode Clean and cold gas 

Electrical power 40 kW 

Gasifier type Downdraft 

Gasification temperature 1050 — 1100°C 

Gas flow 200 m3/h 

Average calorific value 1000 kcal/m3 

Cooling and cleaning system Venturi Scrubber, fat filters, fine filters and safety filter 
Start Through the boot scrubber pump and blower 

Fuel type and size Waste wood with maximum dimensions: diameter 25-50 mm and length 50-60mm 
Fuel consumption 1.3-1.4 kg/kW. Wood containing less than 20% humidity 
Allowed moisture 5 — 20% 

Nominal consumption per hour 60 — 65 kg/h 

Ash discharge 5 — 696 

Normal efficiency conversion >75% 

Normal gas composition N2:50%, CO:1943%, H2: 18+2%, CO2: 102396, CH4:<3% 


The energy generated by the plant (engine-generator) depends on the calorific value and gas flow 
entering and these parameters are a function of the air entering the process in the gasifier [23]. A physical 
analysis of the performance of each of the parts of the system is performed in this paper to determine the 
frequency supplied by the generator. To develop the proposed model of the system, the authors use existing 
models, simulations, and experimental values. 


Table 2. Blower and RICE-Generator Specifications 


Characteristic Value 
Blower Manufacturer United Eng. Works 
Power (HP) 3.0 
Speed (rpm) 3360 
Voltage (V) 215 
Static pressure (mmwc) 700 
Flow (Nm3/h) 300 
Generation Plant Manufacturer Prakash 
Model PN-6 
Power (HP) 65 
Speed (rpm) 1,800 
Alternator Manufacturer Kirloskar 
Power (kVA) 55 
Voltage (V) 220 
Current (A) 145 
Speed (rpm) 1800 
Phases 3 


3. MATERIALS AND METHODS 

The methodology used to perform the global physical modelling of the system and implement the 
frequency control architecture is shown in Figure 3. The first step of the proposed methodology consists in a 
physical modelling of the whole system. The main parts of the system are: the gasifier (the model considers 
static and dynamic behavior), the blower (controls the gas flow), and the RICE-generator set. Secondly, 
validation of the model is performed in open loop using measurements obtained from the real system. 
A cascade control loop is proposed to develop frequency control using the model obtained in the previous 
steps. Finally, results are presented. 
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Figure 3. Methodology 


3.1. Downdraft gasifier physical modelling 

Gasification is a high-temperature thermo-chemical process that converts solid or liquid fuel into a 
gaseous fuel. This technique has been known for centuries, and has been used to obtain a cleaner and 
versatile fuel from unsuitable primary fuels. Gasifiers may be classified into two groups: fluid and fixed bed. 
Two types of fixed bed gasifiers exist according to the direction of movement of the biomass and gases: the 
counter-flow or updraft, and parallel-flow or downdraft. Downdraft is the most widely used technique, 
because it produces a cleaner gas than the other technique, despite having a lower thermodynamic 
efficiency [24], [25]. The physical and chemical processes that a solid fuel particle undergoes during 
gasification are [26]-[29]: 

1. Drying: the process by which solid fuel loses humidity. 

2. Primary pyrolysis: this process is fuel thermal decomposition, releasing gases such as CO, CO2, H20, 
H2, CH4 and many hydrocarbons, some of which have long chains called tars. The carbonaceous solid 
residue is carbon. This process has a complicated energy balance, being slightly endothermic or 
exothermic depending on the nature of the fuel. 

3. Secondary pyrolysis: if temperature conditions permit, the released tars are burnt by an endothermic 
process that releases CO, CO2, H20, H2, CH4 and forms particles of coke. 

4. Reactions of gases: if oxygen is present and the temperature is suitable, there is a combustion of gases 
released during the primary and secondary burnings. In addition, reduction or methanation reactions can 
occur-depending on the temperature. 

5. Coal reactions: residual coal from the above reactions is burnt in the presence of oxygen releasing CO and 
CO2. CO2 reduction to CO and CH4 formation from H2 occurs if the temperature is appropriate. 

Different areas may be defined in a downdraft gasifier in which the above described reactions take 
place as shown in Figure 4. These areas are distinguished by temperature. The quality of the gas produced by 
a downdraft gasifier may be predicted by numerical models and by estimating the chemical reactions using 
the theory of chemical kinetics [28], or with experimental values [23]-[30]. This paper takes into account the 
dynamic behavior of the gasifier. 

Usually in the literature, the experimental values of the studied gasifiers are given in average values 
for periods of minutes for different ranges of stoichiometric combustion or ER (ER is the stoichiometric ratio 
of air needed to produce the complete combustion of wood, where ER-1 this corresponds to 5.22 Nm?/kg of 
wood) [23]. In this work, the gas calorific value (HHV, expressed in KJ/Nm?) for the gasifier studied is 
calculated according to the inflow of air (ma, expressed in Nm?/h) by (1): 


2 ma-1.07 


ma-1.3 .07\2 
HHV(ma) = Nea.) peri sar e 4.8(1 — e~ 6(ma-0.5)) (1) 


Where 0.87<ma<1.89. Figure 5 shows the experimental values and those obtained through (1). Figure 6 
shows the implementation carried out in Matlab-Simulink to obtain HHV. 
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Figure 6. Matlab-Simulink block diagram 


Furthermore, gas flow (mg expressed in Nm3/h) produced by air flow during gasification is shown 


in Figure 7 for the analyzed gasifier [30], [31]. The produced gas flow can then be estimated in a linear 
correlation (2) 


mg(t) = 1.2664ma(t) + 0.137 (2) 
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Figure 7. Gas flow versus air flow for the analyzed gasifier 
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Finally, gas thermal power (P(t) expressed in KW) is given by the product of (1) and (2) and a factor 
of 0.278 to express kJ/h into kW. (3) 


P(t) = 0.278HHV(t)mg(t) (3) 


3.2. Blower model 

Electrical frequency directly depends on the blower, since variations of engine speed directly 
influence the gas flow entering the RICE and the calorific value of the gas supplied by the gasifier. Real tests 
of the dynamic response of a downdraft gasifier have been performed in the laboratories of the Institute for 
Energy Engineering at the UPV. The objective was to model air flow depending on blower speed, by 
measuring the inlet flow with a hot wire flow-meter. The control variable is the speed of the induction engine 
that drives the blower and the air flow at the input of the blower. The transfer function for the airflow value 
following a frequency change in the nominal operating point of the gasifier was obtained experimentally (4). 


0.9616 
s2-0.5236s  .8742 


Gma(s) = [Nm? /hHz] (4) 
The experimental response and the simulated response generated using Matlab-Simulink block for a 
10 Hz step input are shown in Figure 8. 


3.3. RICE-generator physical modeling 

The dynamic behavior of the rotational speed of an alternative internal combustion engine (RICE) 
depends of four variables: the calorific value of the fuel, the air flow, the load to be supplied by the engine 
and the spark advance. The engine may be controlled by changing the flow of fuel, and as shown in the 
previous section, this varies the composition of the gas obtained in the gasifier. To analyze the behavior of a 
RICE when changing the fuel supply, the aspects discussed in the following sections must be considered. The 
engine response at different speeds and different power levels depending on fuel consumption is given by the 
lines of Willans as shown in Figure 9. 

As shown in Figure 9, a linear increase of fuel consumption is produced when increasing the 
supplied power for a given engine speed. Since all lines have the same slope it can be stated that the final 
value of the fuel consumption is directly proportional to the force obtained, therefore: 

1. The engine slows when it receives an increase in load for a given fuel consumption. 
2. Ina steady state of load and speed, it is possible to increase fuel consumption to recover initial velocity. 
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Figure 8. Gas flow obtained for a 10Hz step input in the frequency of the blower 
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Figure 9. Lines of Willans 


3.4. RICE-generator set model 

In the literature there are various models of the dynamic behavior of the speed of a RICE [32]-[36]. 
In [34] propose a simple second-order system that considers the spark advance, the valve openings, and the 
load torque as input variables. In [36] present another model that for its simplicity has been chosen as a 
starting point for this work. To consider the physical phenomena described above, the model has been 
modified to include the valve point, the HHV of the gas and the gas flow as input variables as shown in 
Figure 10. 


Gas Flow [m3/h] 


Valve Point Gas mix(k+1) 
Network Frequency 
Net frecuency (rad/s) 


Network frequency 
Spark Trigger rad/s 


Compression to rpm 


Throttle & Manifold 


HHY [KJ/m3] Gas combustion 


Motor dynamics 


Figure 10. Proposed RICE model 


The accelerator has been modelled as a gas injection which depends of two variables: the valve 
opening point and the pressure at the inlet (5), (6) and (7). 


mq = f (0)g (Pn) (5) 


f(@) = [2.82 — 0.0520 + 0.1030? — 0.00060?][0.025mg — 0.25] (6) 
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where: Mm is the gas mix flow rate into the manifold (g/s); 0 is the valve point (deg); Pm is the manifold 
pressure (bar) and Pau» is the atmospheric pressure (bar). 

The gas pressure at the inlet is modelled as a differential equation in which the difference between 
the gas flow reaching the inlet and the amount of gas entering the cylinder govern the change in the amount 
of gas mass with respect to time. Variation of mass with respect to time in accordance with the law of ideal 
gas 1s proportional to the time derivative of the pressure at the inlet. In addition, the amount of fuel mass 
entering the cylinder is a function of pressure at the inlet and the engine speed. 

As this model is for a four-stroke engine with four cylinders in line, the intake, compression, 
combustion and exhaust occur simultaneously in each of the four cylinders. The torque generated 1s modelled 
from an empirical equation given by [36], which relates the mass of air charged into the cylinder, the gas/air 
ratio, spark advance, engine speed and the HHV of the gas mix (8). 


Tg = [7181.3 + 379.36m,, + 21.91(A/F) — 0.85(A/F)? + 0.260 — 0.00286? + 0.027N — 
0.000107N? + 0.00048No + 2.55e mg, — 0.056?m,,,][0.00025HHV] (8) 


where: Mm is the gas mix flow rate in the cylinder (g); A/F is the air to gas ratio; o is the spark advance 
(degrees before top- dead- center); Tz is the torque produced by the engine (Nm) and N is the engine speed 
(rad/s). 

Once the torque generated by the engine is known, acceleration is obtained by the resistance torque 
and rotational inertia engine by (9). 


J.N 2 Tg - T, (9) 


Where: J is the engine rotational moment of inertia (kg: m?); N is the engine acceleration (rad/s) and Ti is 
the torque produced by the load (Nm). 

The generator is a synchronous machine whose speed of rotation (N) is rigidly linked to the 
frequency of the power generation according to (10). 


n=O (10) 


Where p is the number of pole pairs of the machine [37]. Therefore, voltage frequency is defined by the 
RICE-engine speed. 


3.5. Proposed control strategy 

A frequency cascade control of the blower frequency and the valve point in the manifold has been 
chosen using the global physical modelling of the system described in previous sections as shown in 
Figure 11. Cascade control architecture is used to maximize the benefits by decreasing the energy of the 
control actions and increasing the speed of response [38]. 

The PIDI controller obtains the blower frequency from the speed error of the RICE. The PID2 
obtains the valve point from the blower frequency error. This control acts in a steady state by opening the 
valve. As a result, the blower frequency is constant in steady state and this is desirable because this means a 
constant flow of air into the gasifier and so the gasifier operates at a constant operating point. Constant 
pressure regulation using a mechanical valve ensures that the flow of gas leaving the gasifier is constant 
regardless of the subsequent reduction of flow due the valve point. 
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Figure 11. Control architecture block diagram 


4. RESULTS AND DISCUSSION 

By using Matlab-Simulink was possible to simulate the time response of the modelled system under 
open loop and closed loop control architecture. Because the model is reasonably similar to the real system as 
shown in Figure 12, we have calculated the values of the drivers and tested the model in a closed loop. 
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Figure 12. Model and real system open loop response 


A cascade control was designed to control the frequency of the rural microgrid so that overshoot and 
settling time were minimized. In steady state, the valve point is different for each level of load and the blower 
frequency is set at the nominal point that improves the stability of the gasifier. It can be seen in Figure 13 that 
the model is stable in closed loop, and the network frequency and blower frequency are controlled. The air 
flow into the gasifier is controlled by the valve point. A 10 kW slope on the load requires 120 seconds to 
stabilize completely. 
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Figure 13. Model closed loop response: (a) Network frequency (Hz), (b) Valve point (96), 
(c) Blower frequency (Hz) 


5. CONCLUSION 

Gasification of biomass is projected as one of the best power generation options to supply energy to 
isolated areas in developing countries. Adequate frequency control of the voltage supplied by these systems 
is desirable. An accurate model of the system is necessary to design a reliable control. This paper presents a 
physical model of all the parts of the system: gasifier, blower, RICE and generator, while taking into account 
nonlinearities and the dynamic behavior of the production of gas from a biomass gasifier. A global analysis 
has been performed using previous contributions, experimental data and simulations. 

In this paper the operation of a downdraft gasifier is analyzed in detail and we obtain the transfer 
function that governs operation. Gasification is a process with significant instabilities because of the 
heterogeneity of the internal environment and the heterogeneity of the wood for burning. Calorific value and 
the flow of the gas generated by the gasifier are modelled. Furthermore, a RICE model has been adapted to 
include the gas generated by the gasifier as input. The global model is developed in Matlab-Simulink based 
on a real downdraft gasifier installed in Necocli (Colombia) where the model has been validated by analyzing 
the open loop response. 

Finally, a cascade control is proposed to control the frequency of the rural microgrid, so that 
overshoot and settling time are minimized. With this control architecture, the valve point is different for each 
level of load and the blower frequency is at the nominal point that improves the stability of the gasifier in 
steady state. 
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